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Abstract 

^sO ! The interference of two Bose-Einstein condensates, initially in Fock states, can be described in 

(N . 

terms of their relative phase, treated as a random unknown variable. This phase can be understood, 
either as emerging from the measurements, or preexisting to them; in the latter case, the originating 
states could be phase states with unknown phases, so that an average over all their possible values 

Ctf ■ 

is taken. Both points of view lead to a description of probabilities of results of experiments in terms 
of a phase angle, which plays the role of a classical variable. Nevertheless, in some situations, this 
description is not sufficient: another variable, which we call the "quantum angle", emerges from 

<n : 

t-h , the theory. This article studies various manifestations of the quantum angle. We first introduce 

o 

the quantum angle by expressing two Fock states crossing a beam splitter in terms of phase states, 

00 ■ 

and relate the quantum angle to off-diagonal matrix elements in the phase representation. Then 

G\ : 

we consider an experiment with two beam splitters, where two experimenters make dichotomic 
' measurements with two interferometers and detectors that are far apart; the results lead to viola- 

v : 

tions of the Bell-Clauser-Horne-Shimony-Holt inequality (valid for local-realistic theories, including 
classical descriptions of the phase) . Finally, we discuss an experiment where particles from each of 
two sources are either deviated via a beam splitter to a side collector or proceed to the point of 
interference. For a given interference result, we find "population oscillations" in the distributions 
of the deviated particles, which are entirely controlled by the quantum angle. Various versions of 
population oscillation experiments are discussed, with two or three independent condensates. 
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I. INTRODUCTION 



If two or more Bose-Einstein condensates (BEC) merge, they produce a density interfer- 
ence pattern, as shown by spectacular experiments with alkali atoms [1]. The usual explana- 
tion is that, when spontaneous symmetry breaking (SSB) takes place at the Bose-Einstein 
transition, each condensate acquires a random but well-defined phase. The interference pat- 
tern then exhibits the relative phase. The simplest form of this view involves the use of a 
classical complex variable for each condensate given by 



<lMr)> = vV/?(r)e^ (r) (1) 

where n Q( g(r) are the condensate densities and 4> a ,p{ r ) their phases. Another quantum 
treatment of the problem can be carried out by the use of "phase states," which describe a 
state of two condensates having a known relative phase and a fixed total number of particles 
Q - we will discuss the use of phase states in the next section. For systems containing many 
particles the phase then appears as a macroscopic quantity that has classical properties, but 
takes completely independent random values from one realization of the experiment to the 
next. 

However, Bose-Einstein condensates are naturally described by Fock states, states of 
definite particle number, for which the phase is completely undetermined. Various authors 
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- llO] have shown that repeated quantum measurements of the relative phase of two Fock 
states cause a well-defined value to emerge spontaneously, but with a random value. The 
probability of finding M particles, out of a total of N, at positions rj, • • • Tm (M <C N) is 
shown to be given by 



/7T J\ 1V1 
— II [1 + cos(k • r, + A)] (2) 
-7T Z7T j=l 

where k is the wave number difference between the two condensates. The product in 
the integrand can be interpreted as describing the independent individual measurements 
of position with the interference of two waves of relative phase A, resulting in probability 
[(1 + cos(k • r-j + A)] /2; the A integration expresses that this phase is initially completely un- 
known. Nevertheless, after a series of measurements has been performed (still for M <C N) , 
the product of these probabilities in Eq. is found to peak sharply at some particular value 
A , which becomes better and better defined while the experiments accumulate, but takes 



a completely uncorrelated random value from experiment to experiment. Fig. [T] illustrates 
the peaking effect in the integrand in Eq. (j2J) after 200 measurements (the method by which 
we choose the position values is given in Ref. 
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FIG. 1: The integrand Fjf=i (1 + cos(k • r» + A)) of Eq. ([2]) as a function A after 200 measurements 
at positions r^. These first 200 measurements essentially convert a double Fock state into a state 
resembling a phase state, peaked sharply at some particular value Ao„ which is completely random 
from experiment to experiment. 

Eq. ([2]) is quite capable of describing the interference pattern seen in the MIT experiment 



ll|. Note however that the average over all possible phases makes the phase very similar to 



12j, or to an "element of reality" as defined by 



the integrated variable A in Bell's theorem 

n 

Einstein, Podolsky and Rosen [13] - and we know that this notion combined with locality 
leads to contradictions with quantum mechanics. Eq. can thus be seen as a "classical" 
equation, which is unlikely to be able to describe some truly quantum experiments (for 
instance, it cannot violate Bell's theorem). 

In some conditions, quantum interference effects arise so that the description in terms of 
a classical phase is no longer sufficient; a second angle (or its equivalent) becomes necessary: 
the "quantum angle," which controls the amount of "quantumness" in the results of an 
interference experiment. This article discusses the role of the quantum angle in general. 
While it is possible to carry out such a discussion for the position measurements in free 
space, as in Eq. (j2J), it turns out that interferometers with dichotomic outputs provide 
especially interesting results, for instance in terms of quantum nondocality; this is why 
interferometers will be the central subject of this paper. 

In Sec. [Ill we show how this quantum angle already appears in a very simple situation, 



with one single beam splitter on which two Fock states interfere; we relate the quantum 
angle to phase off-diagonal terms. In Sec. II III we study the effects of the quantum angle in 
an experiment with an interferometer providing dicthotomic results in two different regions 
of space, and leading to violations of the Bell inequalities. But other experiments involving 
directly the quantum angle are also possible. One was suggested to us by the recent article of 
Dunningham et al [if]], who considered the interference pattern of three merging condensates 
and the resulting "phase Schrodinger cat state" formed by the remaining (non-measured) 
particles. In Sec. II VI we consider a simplified version of this experiment with two condensates 
only, which interfere on a beam splitter; among the total of N particles, only M interfere 
and are detected at locations 1 and 2; the remaining are deflected near their sources and 
separately counted in detectors 3 and 4 (m a particles from condensate a, and mp particles 
from condensate For fixed numbers of such particles in detectors 1 and 2, the numbers 
found in detectors 3 and 4, as a function of m a , are found to have an oscillating distribution 
- a "fringe" pattern when plotted over an ensemble of such experiments. We will see that 
this effect, which we call "population oscillations" (PO), involves the interference of two 
peaks in the quantum angle distribution; thus such an experiment would also directly reveal 



the existence of the quantum angle. One can show [17j how these oscillations represent an 
example of quantum interference of macroscopically distinct states (QiMDS), a property of 
quantum mechanics that can verify its validity in large scale systems [lsj . 



II. A SIMPLE INTERFEROMETER 



In the derivation of 1 10 
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|15] , both the classical phase A and the quantum angle A had 
similar origins: conservation rules, which take the form of integrals over these angles. Here 
we show that phase states can also be used to obtain the same results, following a reasoning 
that is similar to that found, for instance, in Ref. [2]. Mathematically, of course, the two 
derivations are equivalent; but, physically, it is interesting to obtain the same results from 
two different points of view. 

We consider the experiment schematized in Fig. [2], where two Fock states with popu- 
lations N a , and Np are emitted by two sources, cross a beam splitter, and interfere in the 
regions of detection 1 and 2. Despite the apparent simplicity of this device we have shown in 
a recent paper that remarkably complex detector distributions can result The double 
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Fock state describing the sources is 



N a ,N ) 



1 



alV |0) 



(3) 




where |0) is the vacuum state and a\ creates particle in state a corresponding to one source 
and ajj creates a /3-state particle corresponding to the other source. The total number of 
particles is N = N a + Np. 



FIG. 2: Two Fock states, with populations N a and Np, pass through a beam splitter, and are then 
made to interfere at detectors 1 and 2. 

The destruction operators a\ and a 2 associated with the output modes can be written in 
terms of the mode operators at the sources a a , dp by tracing back from the detectors to the 
sources, with a phase shift of ix/2 at each reflection: 



The probability amplitude describing the system after crossing the beam splitter with 
mi, 777.2 particles in the detector regions is 



with mi + m 2 = N. To compute the state a™ 1 a™ 2 \N a ,Np), we expand the double Fock 
state in normalized (relative) phase states, defined for two condensates (with constant total 
particle number N) as 






(4) 




(5) 
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(4 + e«4) N |0) 
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E 



N\ 



V2 N N\fo^/n\(N-n)\ 



(6) 



where iV = N a + Nr. The expansion in terms of the phase states is 



The action on phase states of the operators given in (J4j) is particularly simple; if we write 
them as 

en = v ia a ia + Vipap, (8) 
(with Vi a and t> i( g identified by Eqs. (TjE])) we merely obtain 1 Jll|: 



a t |0, N) = J j(v iQ + vve*) \<j>,N-l) 



Applying this result several times to then gives 



where mi + m 2 = A" and 

R(4>) = 2 N ' 2 e- m ^{v la + v w e^(v 2Q + v w e^) 



m 2 a —iN B (j>i 



l + ie 



% + e 



(9) 



(10) 



(11) 



When we insert this result into Eq. (JSj) and take the square modulus, we obtain the 
probability in the form: 



P, 



m\mi 



and find upon multiplying all these factors out: 



n 2ir J_ v 2tt 



(12) 



Pm im2 = 1 ^Py / ^- / ^<x*[(N a -N p ) (0-00/2] 
mi!m 2 ! y_ 7r 2tt 2vr 



x 



+ <j)' -n 



cos 



COS 



/0-0'\ 

cos — COS 



7T 



(13) 



1 The phase state is obtained by repeated actions of the creation operator aj, = (a^+e^Og) over vacuum, but 
no action of the "orthogonal" creation operator aJ, +7r / 2 ' ^ ne ac ti° n of the annihilation operator associated 
to the former operator is therefore simple, while that of the latter gives zero. Expanding the a,'s over the 
a Qji a, and then over and a0 +7r / 2 i and keeping only the component on the first annihilation operator, 
then directly leads to ©. 
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It is then natural to make a variable change by introducing the average of the two phases 



A = (0 + 0' - 7r)/2 



(14) 



now identified as "the phase angle", as well as the difference 



A = (0-0')/2 



(15) 



which we call the "quantum angle." Eq. (fT3|) then becomes 




x [cos (A) + cos (A)] mi [cos (A) - cos (A)] 



(16) 



This probability is a double sum over the variables A and A of a function of these variables 
as well as of the results mi and m 2 . According to (j!5p . if one sets A = or A = ±7T 
in this function, one obtains the contributions of the terms that are diagonal in the phase 
representation. The relevant values of the phase in the initial state then appear directly. 
For instance, if the function has a single narrow peak around some particular value, the 
phase is well-defined; if it has several peaks at various values of the phase, for a pure state 
the system is in a coherent superposition of different values of the phase (a "Schrodinger 
cat" if these values are very different and if the system contains many particles). The role 
of the quantum angle A is precisely to signal the coherent character of the different values 
of the phase (off-diagonal terms in the phase representation). Each time non-zero values 
of this quantum angle play a role, the classical description of Eq. (JSJ) is not sufficient; the 
non-classical behavior occurs because the factors [cos (A) ±cos(A)]/2 in the integrand of 
ffTB"]) can become negative, so that they can no longer be interpreted as probabilities. In the 
(A, A) plane, we will call the "classical region" the region that lies around the A axis at A = 0, 
and the "quantum region" the rest of the plane. 2 

In Fig. EJ^a), we see the absolute square of the coefficient R of Eq. fllljl . showing two 
peaks for a particular choice of mi and m 2 and with N a = Np and M = N. This is not 

2 The regions at A = ±7r are also classical (i.e., equivalent to A = 0), as can be seen by showing that 
the integration segments — tt < A < — tt/2 and tt/2 < A < tt (or equivalently tt/2 < A < 3tt/2) give an 
identical contribution as the region —tt/2 < A < tt/2. To do so make the substitution A' = A — tt and 



A' = 



X-tt. 
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surprising since, classically, an ambiguity in the sign of the phase angle difference also occurs 
in this interferometer: two different values of this difference lead to the same intensities in 
the two output arms. Fig. [3](b) shows a plot of the corresponding integrand of Eq. (flB]) . 
The diagonal phase contributions arise from the peaks on the lines A = 0, ±7r. Here the 
system is in a pure state, so that these peaks are necessarily coherent; peaks in the quantum 
regions (away from A = 0, ±7r.) are also visible, which have a negative sign and therefore 
signal destructive interference (for these particular results of measurement; for other values, 
it is constructive). 




FIG. 3: (Color online) The left side (a) shows the coefficients |-R(0)| 2 in an experiment where one 
finds mi = 9 and mi = 23. The interferometer is not able to distinguish two values of the phase 
difference 4>o between two sources; there remains an ambiguity between ^ + (f>o or ^ — 0o- Here 
00 = ±1.12. The right side (b) shows the integrand of Eq. (|16p for the same values of variables, 
and N a = Ng. Regions on the lines A = 0, ±7r, correspond to diagonal terms in relative phase, 
while those elsewhere represent off-diagonal contributions. 

In Fig. H]we show a particular example of the probability distribution for finding the set 
of {mi, mi} particles in the detectors. The structure has a surprisingly complex dependence 
on the numbers of particles in the Fock state inputs. The simple interferometer is discussed 
more completely in a separate publication |19j . 

In this section we have recovered by the use of phase states the basic results obtained 
from conservation rules in [lol . \\\ . The present method illustrates the relation between 
the two angles and the diagonal or off-diagonal phase terms, and therefore the role of the 
classical and non-classical region in the A, A plane. We now examine how the quantum angle 
changes the description of some other processes for Bose-Einstein condensates involving 
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FIG. 4: The probability distribution of Eq. (Tl6j) for input numbers N a = 26, Na = 24. 



several interferometers. 



III. DOUBLE INTERFEROMETER 



We now discuss the role of the quantum angle in an interferometer experiment d esig ned 
to observe violations of the Bell-Clauser-Horne-Shimony-Holt (BCHSH) inequality [20|, al- 
ready discussed in [14]. The device is shown in Fig. [5] and involves a twin Fock state 
entering a double interferometer, which can be used to measure the relative phase of the 
two condensates in two remote regions of space. The relevance of twin Fock states for phase 

n 

measurements in simple interferometers was already discussed in Ref. |21| in 1993. The 
measurement of the phase of an arbitrary quantum state at different locations of space was 



discussed in Ref. 



22| in 1994. A general discussion of the properties of the quantum operator 



23|. A more 



associated with the phase difference between two modes can be found in Ref. 
recent Re, ,24] gl ves a d.scussioa of the internee of two Foch states and of the detai.s of 
the statistics of the position measurements, in the context of interferences in free space. 



A. Quantum calculation 

For completeness, we briefly recall the quantum calculation in this subsection. The 
destruction operators associated with the output modes can be written in terms of 

the mode operators at the sources a a , ap, a a > and apt by tracing back from the detectors to 
the sources, with a phase shift of n/2 at each reflection, £ or 9 at the shifters, and a l/y/2 at 
each beam splitter. This gives the projections of the two different source modes onto each 
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FIG. 5: Two Fock states, with populations N a and Ng, enter beam splitters, and are then made 
to interfere in two different regions of space Da and Db, with detectors 1 and 2 in the former, 3 
and 4 in the latter. In each of the channels j = 1, 2, 3, 4 particles are counted. We assume that no 
particle is missed: the sum M of the four mj's is equal to N = N a + Ng. 

detector mode 
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ie l< ^a a + tag 
ia a + ie ld ag 
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a a — e ag 



(17) 



where we have eliminated a a > and ag', which contribute only vacuum. The source state, 
having N a and Ng particles in the two condensates is again given by Eq. ([3]). The amplitude 
describing the system crossing all beam splitters with mi ■ ■ ■ 7714 particles in the detectors is 



a 



mi,..,m4 



(mi,m 2 ,m 3 ,m 4 \ $) = (0| 



aT---aT ^ N ^ N " 



a a ag 
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W ■ ■ - m 4 ! yJNJNpl 
The calculation is similar to that of Sec. [TT] and can be found in Refs. 



(18) 



IJ, |l5|. We 



substitute (TTT1) into this expression, make binomial expansions of the sums, evaluate the 
expectation value of the operators, replace Kronecker 5's by integrals in the form 5tv 7 , p = 
e t ( p ~ N ~ , ^ x ~< with 7 = a, /3, and make an appropriate variable change. We then obtain 



V(mi, m 2 , m 3 , m 4 ) 



2- N N a \Ng\ r dx r dA r/Ar AM411 t r . 

i 1 / n~ / — cos (A^q,— iya)A JJ cos A + ^ cos (A — 



i=l 



where r/i = 773 = 1; r/ 2 = r/ 4 = -1; y?i = y? 2 = ~C; ^3 = <^4 = 0. In 



(19) 

15] we also consider 



the case where only M particles are measured among a total number of N, by assuming 
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losses of particles either near the sources or near the detectors. The sum of the probabilities 
associated with the orthogonal states corresponding to the result of measurement is then 

n 2 N ~ 2M M\{N/2\) 2 r dX f w dA r Al v_AfAr a 

V{m 1 ,m 2 ,m 3 ,m 4 ) = — — - — / — / —cos A M cos A + rji cos (A — 

Nl mi! ■■■m A \ J- w 2vr J_ w 2vr fJi 

(20) 

(for simplicity, from now on we assume that N a = N@; the probabilities have now been 
normalized to a total probability of 1 for all events associated with the detection of M 
particles). 

We have associated values of r\ that are +1 for channels of detection 1 and 3, —1 for 
channels detectors 2 and 4. Assume now that Alice, in the first detection region 1, calculates 
the product of all rj values that she obtains, that is the local parity (— l) m2 , which is called 
A = ±1; similarly Bob, in the second detection region 2, calculates B = (— l) m4 = ±1. 
We then have two functions to which the BCHSH theorem can be applied. The quantum 
average of their product is: 

(AB)= (-ir 2+m4 Hmi,m 2 ,m 3 ,m A ) (21) 

The result for the case where all particles are measured [M = N) is found to be fl^ : 

N 

(22) 

B. Classical phase situations 

We consider the case where M N particles are detected; in f)20p . the factor [cos A]^ - *^ 
is peaked sharply at A = 0. Setting cos A to unity in the product and doing the integral over 
A gives 

m\ n dX 4 

V( mi ,m 2 , m 3 , m 4 ) = 4 M mi! .'.. m4 ; J ^ ^ II t 1 + ^ cos ( A ~ ^F* (23) 

where we have taken the — > oo limit of the normalization factor. The quantum angle A 
has now disappeared from the result, so that in the integrand all the terms in the product are 
positive and can be interpreted as probabilities. The BCHSH inequality [20] then provides 

(AB) + (AB') + (A'B) - (A'B'} < 2 (24) 



(AB) 



cos 
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where letters with and without primes imply measurements at differing angles. No violation 
of this inequality is possible as long as (j23|) applies. 

This inequality can also be checked explicitly by computing the value of the average (AB) 
from ( }2"5|) ; we find 



C. Fully quantum situations 

We now assume that all particles are measured. For convenience, Alice's measurement 
angle is taken as (p a = Q/2 and Bob's as 0& = —9/2. We define E(cj) a — 0&) = cos^ (0 a — 0&), 
set 4> a — (fib = (pb — 'Pa' = 4>v — 4>a = £ and cf)b'—(pa' = 3£. We now maximize Q = 3E(£)—E(3£) 
in order to find the greatest violation of the inequality for each N. For N = 2 we find 
Q max = 2.41 at £ = 0.39; for N = A, Q max = 2.36 at £ = 0.26; and for iV -> oo, Q max -> 2.32 
with £ = 0.52/ y/N. The system continues to violate local realism for arbitrarily large 
condensates. 

Despite the identical dependence in the cosine factor in (1251) and ( 1221) . the effect of the 
prefactor, always equal to or less than 1/2 in the classical case, is to prevent the violations 
of the inequalities to occur. Actually, quantum violations disappear even when only one 
particle is missed in the measurement process (M = iV — 1) as discussed in Ref. jljj]. 

D. Discussion 

It is interesting to see in more detail how the quantum angle is involved in the BCHSH 
violation. For instance, Fig. [6] shows the variations as a function of A and A of the function 
that appears in the integral of Eq. ([20]), for N a = N p = M = 40, and 6 = ( = 0. The 
left part of the figure assumes that mi = 6, ni2 = 14, = 14 and WI4 = 6, the right part 
that mi = 6, m 2 = 14, m 3 = 15 and = 5; one immediately notices that, depending on 
the parity of the sum m 2 + m 4 , the peaks in the "quantum region" A 7^ have the opposite 
sign. This explains why the quantum effects will be enhanced if Alice and Bob decide to 
choose the parities (product of all their results 77's) as their local observables A and B. It is 
then natural that strong violations of the BCHSH inequalities should be obtained for this 
particular choice, while of course Alice and Bob could combine their local results in many 




M 



(25) 
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FIG. 6: (Color online) Plot as a function of A and A of the integrand in Eq. (I20p . for N a = Np = 
M = 40, and 9 = £ = 0. Left: m,2 + mi = 20; right rri2 + vn± = 19. Depending on the parity of 
this sum, the peaks in the "quantum region" A / have the opposite sign; this indicates that the 
measurements of parities should be a good choice of local observables to obtain strong violations of 
the BCHSH inequalities 



other ways to obtain functions A and B. 

Suppose now we delete the leading normalization factors in each of Eqs. (120]) and ( )23|) 
and then evaluate the unnormalized values of (AB) for M = 2. The result in each case is 



4 



cos , , 



2 



. Thus the entire difference between quantum and classical averages is in the 
normalization given, respectively, by the integrals over A and A of 

i 

L qu {i, A, A) = t cos A + Vi cos (A + C)P [cos A + rji cos (A + C)P 

mi---m4 

x [cos A + r/i cos (A - 6)} m ' A [cos A + m cos (A - 6)} mA (26) 



where the sum is on all m, totaling 2; and 



A) = [1 + Vi cos (A + C)P [1 + Vi cos (A + C)]™ 2 

x [1 + rji cos (A - #)p [1 + Vi cos (A - 6)} m4 (27) 



For M = 2 we explicitly get 



L gn (£,A,A) = 8 cos 2 A (28) 
Ld(£,\) = 8 (29) 
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The quantum normalization integrand clearly yields a smaller normalization integral en- 
hancing the (AB) average and allowing the violation of the BCHSH inequality. It is this 
variation with quantum angle that allows the violation. 



IV. POPULATION OSCILLATIONS 

Dunningham et al jl6| have considered a situation in which three condensates, a, b, and 
c, each contain initially N/3 particles. A number of them, M < N, form an interference 
pattern on a screen, while the remaining particles m a , m& and m c are counted elsewhere 
(perhaps having been deflected by beam splitters while traveling from the sources), or in a 
second step of the experiment. The numbers of such particles, as a function of m a and mb, 
are found to have an oscillating distribution when plotted over an ensemble of experiments 
corresponding to the same interference pattern for the first M particles. This phenomenon 
was explained as arising from the interference of the two coherent components of a phase 
"Schrodinger cat state" of the system. 



A. Population oscillations by two-source interferometer 

Here we present a simpler version of this experiment based on the interferometer shown in 
Fig. [7J which nevertheless retains the essential features of the three condensate device. The 
general idea is that condensates provide, in a sense, many realizations of the same single 
particle quantum state, since they contain many particles in the same individual state. 
One can then perform experiments where some particles are used to measure one quantum 
observable, some others another "incompatible" observable, which would be impossible with 
one single realization of the quantum state. In this case, the incompatible (non-commuting) 
observables will be the phase and the number of particles. 



1. Experimental setup 

In our version of the experiment, M particles from the two sources interfere in the detector 
D made up of a beam splitter and subdetectors 1 and 2; the other particles are detected 
before they reach the interferometer, with the help of additional beam splitters followed 
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by detectors 3 and 4. For the sake of simplicity, we assume that all beam splitters have 
1/2 reflectivity and transmissivity; we write mi and rri2 as the number of particles seen at 
subdetectors 1 and 2, respectively (with M = m\ +m 2 ), m a and mp the number of particles 
seen in detectors 3 and 4. In this scheme, some of the particles are used to measure the 
relative phase of the two sources N a and Np, the others to obtain information about their 
initial populations. 



FIG. 7: Two source condensates states, with populations N a and Np, emit particles that cross 
beam splitters. Some particles reach the central beam splitter followed by detectors 1 and 2, 
registering mi and m 2 counts. The other particles are then described by a quantum superposition 
of macroscopically distinct states propagating inside the region shown with a dotted line; they 
eventually reach counters 3 and 4, which register m a and mp counts respectively. A phase shift u 
= 7t/2 occurs in one path. 

Assume for a moment that the central beam splitter is removed, so that no interference 
effect between the sources takes place at detectors 1 and 2. Then the experiment separates 
into two independent parts: the detectors 1 and 3 measure the population of one source, 
and the sum m\ +m a gives an exact measurement of the initial population N a ; of course m x 
and m a may fluctuate separately, with a constant sum, but their most likely value is N a /2. 
Similarly, detectors 2 and 4 give information on the population of the other source, and the 
most likely number of their counts is Np/2. 

Now, when the central beam splitter is inserted, the counts of detectors 1 and 2 can no 
longer be ascribed to any of the sources, which are indistinguishable for the detectors; what 
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they actually measure is their relative phase. In classical optics for instance, if the sources 
are lasers with the same intensity and a phase difference <fi a p, the numbers of counts mi 
and m 2 are respectively proportional to sin 2 (A/2) and cos 2 (A/2), with 3 A = (p a p — 7r/2; the 
counting rates therefore provide information about the absolute value A, but not its sign. In 
quantum mechanics, this sign uncertainty becomes an essential ingredient for the creation of 
a superposition of two states with different phases (a "Schrodinger cat"): the measurement 
process at the interferometer projects the initial state of the system onto two categories of 
phase states with opposite phase difference, between which no selection is made. Therefore, 
after the interference measurement, the system reaches a coherent superposition of states 
with opposite values of the phase. 

How can this superposition be observed? The conjugate variable of the relative phase 
is the population difference between the sources; therefore, as the authors of Ref. 16| have 



remarked, if one measures the absolute value of this difference, one expects to see interference 
effects between the two components of the coherent state with opposite signs for the phase. 
Fortunately, even with the central beam splitter inserted, detectors 3 and 4 can still be used 
to obtain information about the populations of the sources. So, for one given value of the 
ratio mi/m2, one expects oscillations of the probabilities associated with given values of 777.3 
and m4, that is "population oscillations". This is the general physical idea, based on the fact 
that Fock states provide many realizations of one single particle quantum state, as mentioned 
in the introduction of this section. We will see that, in our analysis, the interference that 
produces the oscillations occurs between peaks in the quantum-angle distribution. 

Leggett [18] has considered how one might observe coherent superpositions of large num- 
bers of particles by observing their interferenece ("quantum interference of macroscopically 
distinct states" or QIMDS). One can tell the difference between such a pure state and a 
statistical mixture only by observing the off-diagonal matrix elements between the different 
wave function elements. Our population oscillations are the result of such an interference as 
we will discuss below. 

The experimental setup of Fig. [7] is completely defined, as required in the Copenhagen 
view of quantum mechanics; in particular, the setup does not have to be changed from an 
interference setup to a population measurement setup in the middle of the experiment. We 



A phase shift ir/2 is introduced by each reflection on a beam splitters 
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now calculate the probabilities associated with the various possible results of measurements. 



2. Qualitative analysis 

We assume that all particles are detected; the total number then is N = N a + Np = 
mi + ni2 + m a + mp. We will vary the number of particles in detectors 3 and 4 at constant 
N, M, mi,m 2 to examine the behavior of the probability on the set {m a mp}. The destruc- 
tion operators for particles at the detectors in terms of the source destruction operators 
are 

«i = -{a a + iap); a 2 = -{ia a + a p ) 

1 1 

a 3 = ^J^a] a A = ~~^ a P (20) 

The probability amplitude for detecting the set {mi, rri2, m a mp} is given by 

" \N a Np) (31) 

J mi\m2m a \mp\ 

Expand the double Fock state in phase states (Eq. ([7j)) and operate with a™ 1 a™ 2 so the 
state created by the interferometer detectors 1 and 2 is 



\T)^a^ar\N a Np) = ] J^^ jT ge^it^) |0, N - M) (32) 
where m,\ + m 2 = M and 

R^) = (e w + ie^) mi (ie ie + e i(t> ) m2 (33) 
If we take 9 = n/2 then R{4>) takes the simple form 



J4>/2\M 



R(<j>) = (2ie^ 2 ) M cos^ sin| (34) 



Figure [8] shows T((f>) = i?(0)(2ze*^ >/ ' 2 ) _M , which has two peaks at ±0o = ± arctan y 7712/^1 ■ 
This is not surprising since, classically, the ratios of the intensities in the output arms of the 
interferometer determines the absolute value of the phase difference between the two input 
arms but not its sign. Separating negative and positive contributions of provides 

|r> = IV>+> + (-i) m2 IV>-> (35) 
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|T(0)| 




FIG. 8: Variations of R((f>) obtained for mi = 17 and m,2 = 83. The peaks are at 0o = ±0.737r 
(the phase choice 9 = ir/2 gives symmetrical peaks about zero). The relative sign of the two peaks 
is (-l)™ 2 . For large numbers of particles, the measurement produces a coherent superposition of 
macroscopically distinct states ("Schrodinger cat"). 

where 

|^±) ~ e T*(N?-M/2)t | ± O; at _ M) (36) 

We assume that M is large, so that the peaks are sharp and these two branches are orthogonal 
for any 0o not too near zero; and they are macroscopic as long as N — M is large. The 
interference between these two states (QIMDS) is provided by the side detectors in Fig. [7J 
Because we have 

a^a™ 2 | ±0o, iV — M) ~ e ±im ^° |0) (37) 
then the probability of gettting the set {mx,m2,m a ,mp} is 

P(mi, m 2 , m a , mp) ~ 1 + (-l) m2 cos[(m a - m^fio] (38) 

where we have taken N a = Np. The cosine terms in this come from the two cross terms 

. If one does the interferometer experiment for fixed source num- 
bers, say, N a = Np, and considers only those experiments having the same mi, m-i then the 
interference between the two elements will show up in a cosine variation of probability with 
m a . We call this effect "population oscillations." These oscillations are beyond SSB since 
they disappear if one starts from either of Eqs. ([T]) or With a phase state of phase x 
for instance, the action of the destruction operators dip on this state introduces \ instead 
of an integration variable into Eq. (I3"2j) without the integral. No interference effect 
between two phase peaks occurs and the probability is proportional to |i?(x)| 2 . One gets a 
m a ,mp dependence of the probability that is proportional to a simple binomial distribution 
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(N — M)\/m a \mp\, without any oscillation. Actually the angle \ plays no role at all in 
this dependence, which is natural since detectors 3 and 4 do not see an interference effect 
between two beams; they just measure the intensities of two independent sources after a 
beam splitter at their output. 



3. Exact Quantum calculation 



The probability amplitude for detecting the set {m 1; m 2 , m^m^} can be manipulated 
differently: 

1 



n 



^mi\m2m a \mp\N a \Np\ 



NJ.Nr\ 



(0\a? 1 a^c^a? p a]?'*al Np |0> 



E 



mi< 



m 2 \ 



•mi-p-gr r 
^ b 1 u p+q+m a ,N a u rni+rn2—p—q+rnp.Np 



m a +ma 



X 



mi 

E 



2mi+rri2 



(39) 



p=0 P ! ( m l - P)!^ ~ m » ~ P) ! (P + m « + m 2 - ^Va)! 

The probability of getting the set {mim 2 ,m a ,mp} for the sources numbers N a ,Np is 
then 

mi!m 2 !A^o,!A^«! 



P(m 1 ,m 2 ,m a ,mp) 



m a \mp\2 m i +m *2 N 



E 



(40) 



p\(m l - p)\(N a - m a - p)\(p + m a + m 2 - N a )\ 
a result that allows simple numerical computations. 

An alternative form suitable for illustrating the phase relations is obtained if we choose 
to replace one of the ^-functions in Eq. (139)) by an integral, that is 



6 



p+q+m a ,N a 



Zf_J{p+q+m a -N a )4> 

2tt 



The other 5- function simply requires iV = mi + m 2 + m a + nip. For the amplitude we then 

get 

.Y,.LY ( ! [ 



a 



mim2,m a ,m0 



27T 



(41) 
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Squaring C introduces another angle 0'. A change of variables to the relative phase angle 

A = (0 + 0'- tt)/2 (42) 



and the quantum angle 

gives the form 

P(m 1 ,m 2 ,m a ,mp) = 



A = (0-0')/2 



(43) 



NJN l 



w dX r dA 



mi\m2lm a \mp\2 N J_ w 2tt J_ n 2tt 
x [cos A + cos A] mi [cos A - cos Ap 



- i ( N a - m a — +mp ) A 



N a \N p \ 



mi\m2\m a \mp\2 N 
x [cos A + cos A] mi [cos A — cos A] 



r d\ r dA 

/ — / — cos [( N * -m a -N p + mp) A] 

J-7T ^ J-7T 



(44) 



Again we see the appearance of the quantum angle A. We can limit the integration over A 
to non-redundant regions by noting that a segment from tt/2 to 3tt/2 is identical to that 
from just — 7r/2 to tt/2, as seen by making the substitutions A' = A — tt and X' = X — tt. (Cf. 
footnote 2 in Sec. II.) 

A typical example of a population oscillation plot computed from Eq. (14T)]) is shown in 
Fig. |9j We use N a = Np, but from Eqs (14T)]) and (144]) we see that the result would not 
change if we took N a ^ Np. 

Prob 




40 50 60 70 



m„ 



FIG. 9: Plot of P(m 1 ,m a ) given by Eq. (J40]) versus m a for N a = Np = M = 100, mi = 17 and 
m2 = 83. If m2 is even, the central dip is replaced by a peak. 



4- Classical and quantum regions for the distribution 

In Eq. (144]) the m a and mp dependencies appear as a cosine Fourier transform with respect 
to the quantum angle A variable; this cosine Fourier transform is therefore the origin of the 
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population oscillations. If A is set to zero, all m a and mp dependence, and therefore the 
population oscillations, completely disappear. 

We will therefore now concentrate on the distribution F(A, A) that appears in Eq. (PH]) : 

F(A, A) = [cos A + cos A] mi [cos A - cos A] M ~ mi (45) 

and study its variations as a function of the two variables, A and A. As in section [Til the 
band near A = will be called the "classical region", the rest of the A, A plane the "quantum 
region". 

By taking the derivatives of the function F(A, A), we find 4 that the peaks occur at 



A m = and A m = ±2 arctan ( J~~ ) (46) 



A m = ±2arctan ( J— ] and A m = (!, ) 



A m = ±2 arctan ( . / — - ] and X m = ix (48) 

VV m i/ 

The peaks given by ([4*6]) fall in the classical region, and their position depends on the 
observed ratio between m x and m 2 ; this is expected classically since the ratio of the two 
intensities at the interferometer depends on the relative phase of the two inputs. The other 
peaks fall in the quantum region, and will be studied graphically in the next subsection. 



5. Graphical discussion; population oscillations 

We make plots of the quantity .F(A, A) by assuming that N a = Np = M = 40. The mul- 
tiple peaks are visible in Fig. [10] for mi = 17 and rri2 = 23, as well as vii\ = 18 and r«2 = 22. 
The peaks in the figures occur, for m x = 17, at (A, A) = (0, ±1.72), (±1.72, 0), (±1.42, ±7r) 
and, for m x = 18, at (A, A) = (0, ±1.67), (±1.67, 0), (±1.47, ±vr). The two first peaks in the 
"classical region" correspond to Eq. ([4*6]) . while all the other fall in the "quantum region". 
Because peaks corresponding to Eqs. ( 14T]) and f T4"8"]) add up at the border of the diagram, 
for these values of the variables the four quantum peaks at the corners have positions near 

4 A = 7r does not occur because we have eliminated the redundant regions beyond tt/2 < A < 3tt/2 from 
the integral of Eq. £44} . 
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FIG. 10: (Color online) Plots of F(A, A) as a function of A and A for N a = Np = M = 40 with: (a) 
mi = 17 and mi = 23; (b) mi = 18 and mi = 22. Note that — ir < A < it while —ir/2 < A < tt/2. 
In the "classical region" (A ~ ) the two phase peaks have a position that depends on the ratio 
mi/7712, as expected classically. The peaks at the corners of the quantum region change sign with 
the parity of mi and m% they are the source of the populations oscillations shown in Fig. [12] 



7r/2; they are therefore almost independent of the ratio mil mi (if we had chosen smaller mi 
values, these peaks would nevertheless have moved inside 7r/2), in contrast with the classical 
peaks. Moreover, they have a sign that depends on the parity of mi and m 2 , so that it is 
clear that the two kinds of peaks behave rather differently. 

From Eq. (JUJ) we can obtain a quantum angle A distribution given by integrating F: 

— [cos A + cos A] mi [cos A - cos A] M ~ mi (49) 

Here the distribution D(A) has two peaks, as shown in Fig. [HI these peaks are, via the 
Fourier cosine transform, the source of the "population oscillations" as a function of m a . 

Suppose for the moment that the peaks in -D(A) were (^-functions at A = and 7r/2; then 
the cosine transform would be 

/.tt/2 

P(m a ,mp) = / dAcos[(m a -m /3 )A][5(A)±5(A + 7r/2)±5(A-7r/2)] 

J -it/2 

= 1 ± cos[(mQ, — m^)7r/2] (50) 



which oscillates with m a as we have claimed in the form of Eq. ( 1381) . Whether the pattern 
has a maximum or a zero at m a = mp depends on whether m a is odd or even. 



22 




FIG. 11: Plots of -D(A), Eq. (I49p . the integral over A of the function shown in Fig. [TUlfor N a = 

Nfs = M = 40: (a) mi = 17, m 2 = 23 and (b) mi = 18, m 2 = 22. 

The actual plots of P(17, 23, m a , 40 — m a ) and P(18, 22, m a , 40 — m a ) are shown in Fig. 
[T2"j the probability distribution in each case has a finite width, in contrast to the distribution 
shown in Eq. (|50|) . because of the finite width of the peaks in shown in -D(A). The shift 
in phase of the two plots (one vanishing in the middle and the other having a maximum) 
shows that the two components of the interference have changed sign from one case to the 
other. This is precisely the case of the peaks in the quantum region in Fig. [10J Moreover, 
the period of oscillation is constant (maximal for one value of the population and minimal 
for the next), independent of the ratio m\/m 2 , and therefore of the position of the peaks in 
the classical regions. These curves show the results corresponding to the measurements of 
all four quantities m 1) m 2 , • • •, in other words to correlations between various measurements 
at the detectors. If the results are summed over m\ at constant sum m,\ + m 2 , clearly, the 
oscillations wash out. In practice, this means that a post-selection procedure is necessary in 
the experiments. 

When mi < 16 the outer peaks in -D(A) are no longer positioned close to 7r/2 but move 
in to lower A values, and a minimum appears at 7r/2. Nevertheless, oscillations continue to 
occur for values as small as m\ = 1. Only at m\ = does the population oscillation curve 
show just a single central peak. As an example we show the case of m\ = 4 in Fig. [13j 

We finally discuss the A distribution. For this purpose, we sum over variables m a ,m^ to 
get a probability of getting the distribution {mi, m 2 } independent of the source distribution. 
To do the sum we must take into account the relation m Q +m^ = N—M where M = mi+m 2 , 
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FIG. 12: Plots of P(m l7 m a ) of Eq. (gOJ) or Eq. (HP versus m a , for N a 



40 and (a) 



mi = 17, m2 = 23 and (b) mi = 18, ra 2 = 22. Only the integer values of m a are relevant; the linear 
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FIG. 13: Plot of (a) D(A) and (b) P(m 1 ,m a ) for N a = Np = M = 40 and mi = 4,m 2 = 36. 
with M and iV fixed. We obtain (see Appendix A) 

NJNpl r dX r' 2 dA„ (Na _ Ng)A ^ (e iA ) ma (e- lA } 



N-M-m 



P(mi,m 2 ) 



x [cos A + cos A]" 11 [cos A — cos A]™ 2 

NjNpi r d\ r' 2 dA N . M 

— ; ■ »r -, / — / — cos \ (N a — Nr) A (cos A) 



1 7-tt 2vr y_ 7r/2 27r 



m Q !(A — M — m a )! 



tt/2 



x [cos A + cos A]" 11 [cos A — cos A] 



1"<2 



(51) 



The distribution F(X, A) is now multiplied by (cosA)^ - * 7 , which, for large N — M, peaks 
up sharply at A = and damps out all peaks away from A = 0, as shown in Fig. [HJ the 
formula then reduces to the classical case. 
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FIG. 14: (Color online) Plot of (cos A) F(A, A) for N a = Np = 40 and m x = 18, m 2 = 22. Note 
that the peaks of Fig. [Hi] that are away from A = are missing here. 

The classical phase quasi-distribution is then 

— (cos A)^ [cosA + cos A] mi [cosA - cos A] M_mi (52) 

-tt/2 27r 

A plot of this function for the same variable values is shown in Fig. [15j Only the two 
classical peaks survive here. 



Pelass 




FIG. 15: Plot of Pciass(ty of Eq. (|52|) . i.e. the integral over A of the function shown in Fig. [HI for 
N a = Nr = M = 40. The solid line corresponds to m\ = 18, m,2 = 22 with peaks at ±1.67; the 
dotted line is for m\ = 10, rri2 = 30 with peaks at ±2.09. As expected classically, the peaks move 
symmetrically with change in m\. 

An interesting feature of the PO is that, while within the reduced probability of Eq. 
( 15"T|) one can replace A by zero and get the classical limit, it is not correct in Eq. fl4~4"|) . 
which contains no factor (cos A) N ~ M . The result is that one can still get strong populations 
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oscillations and marked even-odd changes, even in the limit M <C N. The quantum angle A 
therefore remains necessary even in this case. 

Population oscillations can continue to exist under certain circumstances even if some 
particles are missed in the measurements; they are more robust in this respect than violations 
of locality. This point is discussed in Appendix A. 



B. Population oscillations with interference fringes in free space 



We now attempt to reproduce the analysis of Dunningham et al. (DBRP) in Ref. [16| in 
which three Fock sources form an interference pattern in free space on a screen, while some 
of the particles are deflected near the sources by beam splitters, where they are counted. Fig. 
[TBI shows the experimental arrangement considered. We will designate M as the number of 
particles involved in the interference measurements made on the screen where interference 
takes place. Then the number of particles measured near the sources having initial particle 
numbers N a = Np = A" 7 = N (as in the work of DBRP) will be m a , mp, and m 7 ; these are 
the particles that did not take part in the interference pattern. All together then we will 
have measured 

3N = M + m a + mp + m 7 (53) 
particles. We can then write the probability as 

P(m a , mp, rrij, r 1; ■ • • , r M ) ~ (T marnprn ^ M T mamfjmj M ) (54) 

where 

a m a a m P Q rn 1 M 

r mam , m7M ) = ° P 7 U( a ^ ka ' n + + a a e ik ^) \N, N, N) (55) 



ijn a \mp\m 7 \ i=i 

To correspond with Ref. 



16] we take k a = k, kp = —k and with k^ = 0. 
We can introduce a vacuum state in between the |r)'s and compute the matrix element 
by multiplying out the interference operators: 

<0|r> = I _ E^ P ^(O(0|a^a^a^a^a^|0) (56) 

Jmalmplm^.y N\ 6 p 

where K p is a coefficient that depends on the r^. The matrix element produces delta functions 
that can be replaced by integrals in our standard way. The results is 

(0 1 r) = I dXa f ^e~ i{ ^~ ma)Aa e~ i{iV ' 3 ~ m ' 3)Aa (57) 

y / m Q ,!m iS !m 7 ! 
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FIG. 16: Particle beams from three sources emitting N a , Np, and iV 7 particles, respectively, interfere 
in free space and can produce an interference pattern on a screen within the central region. Some 
particles (m a , mp, and m 7 ) are deflected near the sources to be counted in detectors. In our 
simulation we have N a = Np = iV 7 = N; N particles reach the interference region and m a + mp + 
7717 = 2N particles are deflected near the sources. 



M 



n(e<* 



1) 



If we take the absolute square of this we introduce two new variables X' a and A^ . We 
then make the following variable changes: 



A, 



A n 



->4 

w A' x A' 



(58) 



The probability then becomes 
P(m a ,mp,m 1 ,r 1 , ■ ■ -,r M ) 



N\ 3 



n dx r dx' r dA r <im 



m a \mp\m^\ J _ w 2n J_ n 2n J _ n 2n J_^ 2tt 



■i(N— m a )A 



xe 



M 

-i(N-m p )A' -Q ^ + e iA + e iM 
i=l 



+2cos(2A;-r J - A-A> i(A+A ' )/2 

+ 2 cos(A; • r, - X)e lk ' 2 + 2 cos(A; • n - X')e tA ' /2 



(59) 



If we sum out the m 8 we find the probability for the r— set under arbitrary source number 
detections: 



P(r 1 ,---,r M ) 



N\ 3 
(37V - M) 



r dx r r d\ r dvv 

! /_„ 27 /_„ 2^ /_„ 2^ /_„ 27' 



-iiV(A+A) 
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,\ 3JV-M 



X 



X 



(l + e iA + e if 

M 

]J [l + e tA + e lA ' + 2 cos(2£; • r< - A - A')e t(A+A ' )/2 



i=l 



+ 2 cos(£; ■ n - X)e iA/2 + 2 cos(fc • n - X')e iA ' 



/2' 



(60) 



The integral method above is not very useful for simulations. We have developed a 
recurrence method in which the wave function for R measurements is written in terms of 
that for R — 1 measurements. We do not give details here to save space. All the probabilities 
p(m a , mp) for finding particles in the source detectors for a given set of positions rq, ■ ■ • , tm 
(with m 7 given by Eq. (1531) ) are computed in a single recurrence run. We show a plot in Fig. 
[T71 of the resulting population oscillation ridges. DBRP found that the ridges were parallel 
to one axis shown in their Fig. 5. Indeed our ridges are parallel to m 7 =constant, which 
would have been more obvious had we plotted using, say, the m a , m 7 axes. The parallel axis 
in our case is the one having the intermediate vector (fc 7 = 0) so that we are in agreement 
with the results of DBRP. 
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FIG. 17: (Color online) The probability of measurement of populations, made with beam splitters 
near the three sources for M = 100 position measurement in the interference pattern in free space 
for three-condensates, where initially each source had N = 100 particles The two horizontal axes 
are the source populations m a and mp of the condensates with opposite wave vectors; the source 
count m 7 , corresponding to the condensate with zero wave vector, is equal to 50 — m a — mp. 
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V. CONCLUSION 

In many cases, such as for instance the description of the MIT experiment [1] with 
initial Fock states, introducing a classical relative phase angle A is sufficient. Cases exist, 
nevertheless, where the classical phase is not able to explain all quantum predictions, and 
where introduction of the quantum angle A (or its equivalent) becomes necessary. We have 
discussed two examples, in Sees. II III and IIVt where interesting physical effects can not be 
understood only in terms of the classical phase. In both cases, quantum effects are related 
to peaks of the function F(A,A) in the "quantum region" (i. e., away from A = 0), and 
disappear completely if A is set to zero. 

In the first double interferometer experiment, we find violations of the BCHSH inequalities 
and therefore violations of locality. Setting the quantum angle to zero reduces the equations 
to purely classical equations, which could be interpreted as being integrated over a hidden 
variable as in Bells theorem. Only the quantum angle leads to the violations. 

In the population oscillation experiment, we find that simultaneous measurements of 
"non-commuting variables" phase and particle number within the same apparatus yield os- 
cillations in measurements of the number variable that are a direct result of the off-diagonal 
phase (i. e., quantum) peaks and provide an example of QIMDS. However, as discussed in 
Appendix B, one can replace the measurement of the phase by that of the parity, which 
does not fix the relative phase of the two condensates at all; but this does not completely 
cancel the population oscillations since the central dark fringe remains present with a 100% 
contrast, while the characteristics of the superposition are completely changed (the "phase 
cat" becomes completely "blurred"). The fact that some population oscillations remain visi- 
ble, at least for the first fringe, illustrates that the PO can exist more generally than with 
just the coherent superpositions of different macroscopic phases. 

The two experiments we have discussed are of somewhat different nature. The former 
exhibits strong quantum non-locality effects, while for the latter we have not found violations 
of the Bell inequalities. Nevertheless, while for the former the violations of the inequalities 
require that all particles are measured (they disappear as soon as a single particle is missed), 
the population oscillations are a manifestation of the quantum angle that is more robust as 
we show in Appendix A; they can still exist, although in a more limited way, when a few 
particles are missed. 
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APPENDICES 

A. Incomplete measurements in the PO experiment 

1. No phase measurements 

We study the experiment of Fig. [7] again, but now assume that no measurement is per- 
formed in the interference region D, and that only the population measurements are per- 
formed; then, whether or not a beam splitter is used in this regions does not matter any- 
more). We then have to sum the probabilities (144)) over m x and m 2 , with a constant sum 
mi + rri2 = M. The summation introduces the M-th power of a binomial [e tA + e~ tA ], but 
only one term of this power survives the A integration; we then obtain 

N 'A«> M> 

m ^M P{mi ' m2 ' m " m " ) = mJr^pKN-py. (61) 

with p defined by 

2p = N a — m a — Np + mp (62) 

(one can easily check that the right hand side of this equation is an even number). The 
A integral has now disappeared, as expected since no measurement of the relative phase is 
made. Moreoever, the probability factorizes as expected since, in the absence of interference 
measurements, two completely independent experiments are performed in different regions 
of space: in each region, the transmission or reflection of the particles on the beam splitter 
are independent random processes. 
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2. No population measurements 



Conversely, assume that all population measurements are ignored and that only the 
interference measurements are considered. The corresponding probability is then 

p, > NjN f \ r d\ r dA K _ M 

m „ + ,£ N -M P(mi ' m2 ' m - ms) = MwLi/-> |mA1 

x [cos A + cos A] mi [cos A -cos A] m2 (63) 

Now the phase A no longer disappears, but combines its effects with the quantum angle A; 
the [cos A] N ~ M introduces a peaking function around the origin, which may behave similarly 
to a delta function if N — M is sufficiently large. We now discuss the interplay between the 
classical phase and the quantum angle A. 



3. Missed particles 

Next suppose some of the particles are lost and not measured in either interferomenter nor 
side detectors of Fig. [7J We have seen in the case of the double interferometer Bell- violation 
experiment that a single missed particle can remove any locality violations. We simulate 
these lost particles in the PO experiment by putting additional side detectors as shown in 
Fig. [THJ Assume that the beam splitters at detectors 5 and 6 each have a transmission 
coefficient T. 

We assume that particle losses m 5 and m§ are known to total Ml, but the individual 
numbers are not actually recorded. Thus to get the probability we are interested in we must 
sum over all and mg adding to the total Ml- Proceeding as in Sec. IV we find 

njnb\t n - Ml (r) Ml r dx r dA r/Ar 

P(m 1 ,m 2 ,m aj m l3 ) = — - - - , ojV _ 2Mr , / 7T / TT cos ~ m ° ~ N P + m P) A 

x cos(A) Mi [cos A + cos A] mi [cos A - cosAp (64) 

with R — 1 — T. The result of the lost particles is the factor cos(A) Mi , which, if Ml is 
large enough, diminishes the quantum peaks, as we have seen before as in, say, Eq. (l63p . 
The result maintains the same form if we also allow particles to be missed elsewhere in the 
device, say, after the beam splitter at detectors 1 and 2. 

If we count Nd = mi + m 2 + m a + mp particles in the real detectors but Ml were missed 
in one place or another in the device, then we must have had = No + Ml particles 
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FIG. 18: Modified population device to show the effect of particles missed in the measurement. The 
missed particles are supposed to enter detectors 5 and 6. We assume we know that + mg = Ml 
particles are missed, but we do not know their distribution in the new side detectors. Thus we sum 
over all m^,m% to get the resulting probability when some particles are not detected. 

in the sources originally. The missed particles could have come from source a or source 
P. We assume that the sources originally have N a = N a + A a and Np = Np + Ap where 
N a + Np = No and Ml = A a + A^. We first fix Ml and sum over all possible A a and 
then sum over all Ml in principle from to oo. If R is small, then the sum over Ml should 
converge after a reasonably small number of missed particles. That is, the probability of 
missing Xl particles, where Xl is very large is negligible. One would hope that if R is 
small enough, then the PO will converge to a situation in which the fringes are not lost. We 
find this to be the case under certain conditions. We can also find the average number of 
particles lost by multiplying the probability by Ml and summing over all A a , Ml, and m a 

Consider the situation with m\ = 17 and rri2 = 83. The PO for the cases with T = 0.98 
and T = 0.99 are shown in Fig. [19j For T = 0.98 the oscillations are completely removed 
and for 0.99 only a remnant is left. In the later case we have lost 3.8 particles on average. 

The smaller the value of m 1; the closer in to A = are the off-diagonal peaks in F(A, A) 
so that they get less blotted out by the cosA Mi factor. If we lower mi to 3 we find the 
results in Fig. [201 We still lose the central dip in the PO diagram for T = 0.98, but for 
T = 0.9 9, we get a much deeper remnant. 

To what degree must we restrict losses to guarantee that we would have more than a 
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FIG. 19: The PO when the total number of detected particles in 200 with m-i = 17, = 83, 
and T = 0.98 (left) and 0.99 (right). Note the very small depression in the center of the right 
plot showing a remnant of the PO after all possible losses are considered. The average numbes of 
particles lost here are 6.5 and 3.8, respectively. 
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FIG. 20: The PO when the total number of detected particles in 200 with m\ = 3, mi = 97, and 
T = 0.98 (left) and 0.99 (right). A much larger remnant depression remains in the 0.99 case here. 
The average numbers of particles lost here are again 6.5 and 3.8, respectively. 



single dip? Fig. |2~T1 shows the case of N = 200, mi = 17 with the transmission coefficient up 
to T = 0.997. Only 1.4 particles have been lost here. For smaller mi values one gets deeper 
central dips, but not the dips on the side for the same T values. 



B. Measuring the parity 

Fig. [10] shows that the A peaks in the quantum region have a sign that depends on 
the parity of mi or m^-, which suggests that a possible method to observe the population 
oscillations is to associate them with a measurement of the parity at the interferometer, 
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FIG. 21: The PO when the total number of detected particles in 200 with mi = 17, m?, = 195, and 
T = 0.997 . Here we get more than just the central dip. The average numbers of particles lost here 
is just 1.4. 

instead of the relative phase of the two condensates. 

Fig. [22] illustrates what is obtained if, for instance, one adds the probabilities associated 
with all odd values of m\. The left part of the figure shows the variations of F(A, A), the 
right part the associated population oscillation as a function of mi, with m-i = 40 — mi. One 
notices the disappearance of the two peaks that characterized the coherent superposition of 
two values of the relative phase; they are now replaced by a more delocalized structure, 
similar to a ridge. In other words, the "Schrodinger cat" is now spread over many values of 
the phase. But one also sees in the right part that the populations oscillations still exist, 
with a central dark fringe that has 100% contrast when m a varies only by one unit; the 
variation is actually not very different from the right part of Fig. [131 except of course the 
change of sign due to the change of parity of mi. This shows that the central fringe of the 
population oscillations is not specifically related to a measurement of the phase, or to the 
existence of any "Schrodinger cat"; it continues to exist if a very different physical quantity is 
measured, such as the parity, which does not give any particular information on the relative 
phase of the two condensates. 

C. Two phase measurements 

In the population oscillation experiment, we have considered the use of a single interfer- 
ometer; we can generalize this experiment to two interferometers each of which have different 
settings, ( and 9. We begin with the device in Fig. [5] and add to that the side detectors 
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FIG. 22: (Color online) Left: plot of F(A, A) as a function of A and A for N a = Np = M = 20, 
obtained by summing the probabilities of all odd values of m-i and mi- One notices that the A 
peaks in the classical region are now spread over many values of A, except sharp variations around 
A = and A = ±7r; the function still takes significant values in the quantum region A ^ 0. Right: 
corresponding population oscillations; the narrow central fringe is still perfectly visible with a 100% 
variation when m\ varies by only one unit only (the central fringe is dark because m\ is odd). 



shown in Fig. [7J to allow phase-type measurements in two different regions of space as well 
as population measurements near the two sources. The resulting apparatus is shown in Fig. 



The calculations proceed as in previous sections and lead to a result for the probability 
of finding the series of results {mi,m2, m 3 , m^, m a , mp} equal to 

NJNpl 



V(mx, m 2 , m 3 , m 4 , m a , mp) 



2 M N\m 1 \---m A \m a \mp\ J_ n 2tt J n 2tt 

x cos [(N a — Np — m a + mp)A] 



x 



J[ [cos A + r]i cos (A + (pi)] mi 



(65) 



i=l 



where M = m 1 + h m 4 , 771 = i] 3 = 1; rj 2 = rj 4 = -1; ip x = cp 2 = —(] cp 3 = ip 4 = 9. This 

result is essentially the same as Eq. (12TJ|) with the addition now of the cosine transform in 
A. This factor allows population oscillations as in Sec. IV. However now we have the option 
of adjusting relative phases between the two interferometer sets. 

A summation version of the probability is much more convenient for computing popula- 
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FIG. 23: Two Fock states, with populations N a and Np, enter beam splitters, and are then made 
to interfere in two different regions of space Da and Db, with detectors 1 and 2 in the former, 3 
and 4 in the latter. In each of the channels j = 1,2,3,4 particles are counted. The extra detector 
5 and 6 count particles that are not measured by the interferometers. 

tion oscillations. A result analogous to Eq. (|4"0]) of Sec. [TV] is 

mj ! . . . m J (—I )P2+P4 e -«(C+<?)(P3+P4) 



V(m 1 ,m 2 , m 3 , m 4 , m a , m p ) 



in !,,, !->.V + 2.\/ E 

1 



ftjWP^ 7 ^ -P2)K m 3 - P3)K m 4 - PaV- 



(N a -m a -p 2 -p3~ Pa)\ 
1 



(66) 



(mi - iV a + m a + p 2 + p 3 + P4.Y 

Because two independent settings 9 and ( are now available, the phase sign ambiguity 
can be removed. As a consequence, by adjusting the phase angles on the interferometers, 
we can now control the relative sizes of the two classical peaks, i. e., those along A = 0. 
Consider the following plots where we show the last line of the integrand of Eq. (IrJo]) and 
the population oscillations given by Eq. (166]) associated with the same parameters. With 
the phase shifters set at zero the two classical peaks have equal sizes and there is a definite 
population oscillation structure (Fig. |24"|) . However, with a different phase shift, one of the 
classical peaks can be made much smaller as seen in Fig. [25] and the quantum peaks become 
smaller as well. Moreover, the population oscillation central zero no longer vanishes. For 
other phase shift angles (for instance C = and 9 = 2.5 radian), the integrand can be reduced 
to a single classical peak with no peak in the quantum region at all; the corresponding 
population oscillation central dip then becomes a simple peak. 

If the state vector is the sum of two components centered around two different values of 
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FIG. 24: (Color online) Double interferometer results for {mi, mz, m^, m^} = {2,8,1,9} with phase 
shifts £ = = 0. Left: The last line of the integrand of Eq. (|65p . There are two equal classical peaks 
(on A = axis) and negative quantum peaks. Right. The corresponding population oscillations 
with a "dark fringe" at the center. 




FIG. 25: (Color online) Double interferometer results for {mi, mz, m^, m^} = {2.8,1,9} with phase 
shifts C = 0, = 2.9. Left: One of the classical peaks has almost completely disappeared and 
the negative quantum peaks are now very much smaller. Right: The central depression of the 
population oscillation plot no longer goes to zero but becomes an indentation. 

the phase, and if the norm of one component is larger than that of the other, one obtains two 
peaks in the classical region, one large and one small. The small classical peak corresponds 
to a population in the phase representation, so that it is second order with respect to the 
second component of the state vector. By contrast, the peaks in the quantum region are 
first order, since they correspond to off-diagonal matrix elements. As a consequence, when 
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one reduces the small phase component, the small classical peak disappears more rapidly 
than the quantum peaks. This explains why the left of Fig. [25] has a classical peak that is 
barely visible, but still clearly shows the (negative) quantum peaks. 
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